Shape-induced phenomena in the finite-size antiferromagnets 
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It is of common knowledge that the direction of easy axis in the finite-size ferromagnetic sample is 
controlled by its shape. In the present paper we show that a similar phenomenon should be observed 
in the compensated antiferromagnets with strong magnetoelastic coupling. Destressing energy which 
originates from the long-range magnetoelastic forces is analogous to demagnetization energy in 
ferromagnetic materials and is responsible for the formation of equilibrium domain structure and 
anisotropy of macroscopic magnetic properties. In particular, crystal shape may be a source of 
additional uniaxial magnetic anisotropy which removes degeneracy of antiferromagnetic vector or 
artificial 4th order anisotropy in the case of a square cross-section sample. In a special case of 
antiferromagnetic nanopillars shape-induced anisotropy can be substantially enhanced due to lattice 
mismatch with the substrate. These effects can be detected by the magnetic rotational torque and 
antiferromagnetic resonance measurements. 
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I. INTRODUCTION 



The fact that antiferromagnetic crystals break up into the regions with different orientation of antiferromagnetic 
vectors below the Neel temperature was predicted theoretically by L. Neeli and then proved experimentally (see, 
e.g. 4——— and many others). 

Domain structures observed in different antiferromagnets have some common features that we summarize below. 

Magnetic domains with different orientations of antiferromagnetic vector are characterized by different tensors 
of spontaneous strain and so can be treated as deformation twins. 

The morphology of antiferromagnetic domains is similar to the morphology of deformation twins in martensites. 
In contrast to ferromagnets, domain structure in antiferromagnets is regular, periodic and consists of alternating 
stripes with different deformation. 

Unlike ferromagnets, domain walls separating domains with nonparallel antiferromagnetic vectors are plane-like 
and are parallel to low-index atomic planes. 

Deformation does not map the orientation of antiferromagnetic vector locally (e.g., inside the domain wall orien- 
tation of antiferromagnetic vector is determined by competition between exchange interaction and deformation- 
induced anisotropy). 

Antiferromagnetic domains spontaneously appear below the Neel temperature. The domain patterns observed 
during heating-cooling cycles through the Neel point may be either identical or similar to each other. 

Domain structure may be reversibly changed by external magnetic field or stress. 

The properties (fi|)- (|iv]) show that magnetoelastic coupling plays the leading role in formation of the domain structure 
in antiferromagnets. It follows from (jvj), (|vT]) that the domain structure may be considered as thermodynamically 
equilibrium notwithstanding the fact that formation of the domain walls is associated with positive contribution into 
free energy of the whole sample. Regularity of the domain structure (properties (JTiJ) , (JTTTJ) ) excludes an entropy of domain 
disorder as a factor which leads to a decrease of free energy of a sample and favors formation of inhomogeneous stated. 
Properties (|Tv|) - (|vT)) may be explained by presence of the elastic defects (dislocations, disclinations, etc.)^ that produce 
inhomogeneous stress field in a sample. This "frozen-in" extraneous (with respect to an ideal crystal) field stabilizes 
inhomogeneous distribution of antiferromagnetic vector via magnetoelastic interactions and ensures reconstruction of 
the domain structure during heating/cooling cycles. 

Another models consistent with all the above mentioned properties is based on the assumption that antiferro- 
magnetic ordering is accompanied by appearance of so called quasiplastic stresses 11 coupled with the orientation of 
antiferromagnetic vector. We assume that these intrinsic stresses are caused by virtual forces that represent the change 
of free energy of the system with displacement of an atom bearing magnetic moment. Self-consistent distribution of 
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the internal stress field depends upon the shape of the sample and is generally inhomogeneous. Equilibrium distribu- 
tion of antiferromagnetic vector maps the stress field and thus is also inhomogeneous and sensitive to application of 
external field and temperature variation. 

Both the defect-based and defectless models exploiting magnetoelastic mechanism predict similar dependence of 
macroscopic characteristics of a sample vs external magnetic and stress field, but lead to different results when applied 
to a set of different samples. Namely, in the framework of the defect-based model, the domain distribution, domain 
size, and some other quantitative characteristics may vary depending on technological conditions and prehistory of 
a sample. On the contrary, the defectless model predicts variation of macroscopic properties of antiferromagnetic 
crystals with variation of their shape. 

Below we predict shape-related phenomena in antiferromagnets that can be experimentally tested. In the frame- 
work of the defectless model we calculate the effective shape-induced anisotropy which can be determined by torque 
measurements, and the frequency of the lowest spin-wave branch detectable by antiferromagnetic resonance technique. 
We consider the case of an "easy-plane" antiferromagnet typical example of which is given by NiO, C0CI2, or KC0F3. 

II. DESTRESSING ENERGY 

According to our main assumption, an antiferromagnetic vector L(r) may be treated as a quasidefect that produces 
intrinsic stress field a^ ln \r). Thus, the thermodynamic potential of the finite-size antiferromagnet can be represented 
as (see, e.g., Ref. [To|) 

$[L(r)] = J |/ raag [L(r)] - i<j (in) [L(r)] : c' 1 : a (in) [L(r)]| dr + $ dcst . (1) 

Here / ma s is a "bare" magnetic energy, the second term is a self-energy of quasidefect, cr 1 is 4-th rank tensor 
of elastic stiffness, <£> dest is a destressing energy^ which describes interaction between the quasidefects localized at 
different points. Sign ":" is used to denote an inner product between the 2nd rank tensors. 

Explicit expression for destressing energy is obtained from the requirement for mechanical equilibrium with due 
account of boundary conditions at the sample surface. The main nonnegative contribution arises from averaged (over 
the sample volume V) internal stress (<r m ) and can be represented as 

$ dc S t = V {a ^. kim (c7 i„ m)) a = j ^ (r _ r , } ^ (2) 

Z OTkOT m Jy 

where G^ rn {Y — r') is a 3D Green's function of elasticity (with zero nonsingular part) and 4-th rank symmetrical 

destressing tensor H depends upon the sample shape. 

Functional dependence between intrinsic stress tensor and antiferromagnetic vector is given by a constitutive relation 
which should satisfy the principles of locality, material objectivity and material symmetry. The simplest form of such 
a relation which assumes isotropy of magnetoelastic properties of the media is 

0% = ^L 2 6 jk + A' (ljLu ^-8 jk )j , (3) 

where coefficients X v and A' define the principal stresses of magnetoelastic nature. 

Substituting © and © into (jTJ) one comes to a closed expression for thermodynamic potential, minimization of 
which gives equilibrium distribution of L throughout the sample. 

Two terms of magnetoelastic origin in JT]) have one principal distinction. The structure of the local energy contri- 
bution (second term) is defined by crystal symmetry, while the structure of $ dost depends upon the sample shape. 
In the framework of phenomenological approach local energy contributes to effective anisotropy constant only, while 
destressing energy is responsible for the domain structure formation and may be a source of artificial 4th order 
anisotropy as will be shown below. 

III. APPLICATION TO AN "EASY-PLANE" ANTIFERROMAGNET 

To understand the role of destressing energy in the shape-induced phenomena we consider the simplest case of an 
"easy-plane" antiferromagnet (point symmetry group of the crystal includes 3-rd, 4-th or 6-th order rotations around 
Z-axis) cut in a form of an elliptic cylinder with a and b semiaxes (parallel to X and Y axes, respectively) and 
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Table I: Shape-induced _K" ; e s las and magnetic f^^,i aae anisotropy (in erg/cm 3 ), and critical aspect ratio for typical "easy-plane" 
antiferromagnets (details of calculation see in^). 

Crystal /™ g planc (a/b) a 

NiO 0.8-10 4 288 1.1 

C0CI2 5.6-10 5 < 3 ■ 10 4 3.4 

KC0F3 3-10 6 5-10 5 1.5 



generatrix parallel to Z. The elastic properties of the media are supposed to be isotropic (en — C12 = 2C44). In this 
case nontrivial contribution to destressing energy takes a form 

$ dcst = ^{Kf"(L$-L%)+K£ m [(I% c -L 2 Y ) 2 

+ 4(L X L Y ) 2 } - Kf™[{L\ Ly) 2 ML X L Y ) 2 }} , (4) 
where effective shape-induced anisotropy constants are 

, ola8 _ a -fe(A') 2 (2-3^) + A„A' 



K c 2 



a + b 4c44(l — v) 



Kf . = (ffl^M, s£= y' « , (5) 



8044(1-^) ' \a + bj 6c 44 (l 

and v = Ci2/(cn + C12) is the Poisson ratio. 

Magnetic energy density of such an antiferromagnet in the external magnetic field H (low compared with spin-flip 
value) may be written as: 

/ mas = \kT & L% + /r S P ia„c - \x[K x L] 2 , (6) 

where x is magnetic susceptibility, out-of-plane anisotropy constant K™ ag ^s> K^ &s is large enough to keep antiferro- 
magnetic vector in XY plane, and explicit form of the in-plane magnetic anisotropy /j^pi^e is specified by a crystal 
symmetry. 

For a typical "easy-plane" antiferromagnet /™ a pi anc is much less than the effective constants (J5J) of magnetoelastic 
nature (see table Up. So, in such a crystal destressing effects may stimulate formation of the domain structure and 
change equilibrium orientation of antiferromagnetic vector. 

It should also be stressed that three different shape-induced anisotropy constants ([5]) depend on the aspect ratio 
a/6 in a different way. This opens a possibility to control macroscopic properties of the sample varying its shape. 



IV. FORMATION OF THE EQUILIBRIUM DOMAIN STRUCTURE 

If in-plane magnetic anisotropy is small but not vanishing, then, the destressing energy favors formation of the 
domain structure. For example, in the case of an isotropic sample (a = b) the only nontrivial term with if. elas in 
eq. is nonnegative. In the absence of external field it can only be diminished by zeroing average values of (L 2 — L 2 ) , 
(L x L y ), i.e., by appearance of equiprobable distribution of domains with different orientation of antiferromagnetic 
vector. 

The external magnetic field causes rotation of antiferromagnetic vector and removes degeneracy of various domains. 
Due to a long-range character of elastic forces, the field-induced ponderomotive force that acts on the domain wall is 
compensated by the destressing, restoring force. Competition of these two factors determines equilibrium proportion 
of different domains. If, for example, external field is applied in-parallel to an "easy" direction in XY plane (say, X 
axis), then, the volume fraction £ of energetically preferable Y-type domain (in which L _L H), increases quadratically 

with field £ = 0.5[1 + (H/Hmd) 2 ] up to monodomainization field Hmd = \J K^ as ' /\ oc A' / ' ^/x c -a- 

The domain structure may be also observed in samples with small but nonzero eccentricity (a « b) providing that 
in-plane magnetic anisotropy is large enough to keep two different equilibrium (stable and metastable) L orientations: 
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Kf 3 * < /™^p lano - In this case it is the shape factor that removes degeneracy and hence equiprobability of the domains. 
In the above example the fraction of X-type domain (L||X) depends on the aspect ratio as follows 



1 Kf as 



2 ^4°an 



and the domain structure reproduces the orthorhombic symmetry of the sample. 

Critical values of aspect ratio (a/b) CI (obtained from the condition if| las = /^pi^J at which equilibrium domain 
structure is still thermodynamically favorable are given in the last column of table Hi 



V. TORQUE EFFECT 

If the aspect ratio a/b of the sample noticeably differs from 1, then, all the effective anisotropy constants in (j4|) 
have the same order of value and are much greater than in-plane anisotropy, S> /^pi ano (see table J}. So, the 

sample has shape-induced uniaxial anisotropy regardless of its crystallographic symmetry. 

An appropriate tool for measuring anisotropy constants is a rotational torque of untwinned crystal in the magnetic 
field. If the rotational axis is perpendicular to "easy-plane" XY and magnetic field makes a ip angle with X axis (see 
inset in fig. [!}, then the rotational torque can be calculated as —d^/dip where free energy potential $ is given by 
eq. |T]). With account of the eqs. (H])-© the rotational torque per unit volume is represented as 



T(ip) = #f las sin 20ft/; 



t>) + 



2Kf™sm49^), 



(8) 



where a 9 angle between antiferromagnetic vector and X axis unambiguously determines equilibrium orientation of L 
and is calculated from the condition for minimum of the potential {1} : 



Kf as sin 26» + 2Kf™ sin 40 - - X H 2 sin 2(0 - V) = . 



(9) 



Analysis of eqs. 



shows that i) effective anisotropy is determined by shape-dependent (via a/b ratio) constants 
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Figure 1: Rotational torque of untwinned NiO crystal at RT in H=4.8 kOe for [111] rotational axis. Triangles - experimental 



data" 
effect). 



solid line - theoretical approximation according to (p|. ((9|, dashed line - approximation of infinite crystal (no shape 



Kf as 



JQan ! u ) shape- induced anisotropy removes multiaxial degeneracy of equilibrium orientation of antiferromag- 
netic vector and thus excludes formation of domain structure; iii) in the field absence the preferred orientation of 
antiferromagnetic vector coincides with the longer ellipse axis (i^2 las > 0, a > 6); iv) external magnetic field applied 
along the longer ellipse axis may induce spin- flop transition at H = Hsf, where spin-flop field 



H S F = 



K c 2 



las 



4an 



A' 



VXC44 V a + b 



(10) 



is governed by the sample shape and is independent of crystalline anisotropy. 
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Figure 2: Rotational torque (calculated) of untwinned NiO crystal at RT in H=5 kOe for [111] rotational axis for the samples 
with different shape. For ellipse-section samples an aspect ratio is a/b = 20 (a); 3 (b); 1.2 (c); (d) - square-section sample. 



An interesting result is obtained for the sample having a square cross-section. In this case equilibrium distribution 
of antiferromagnetic vector is in principle inhomogeneous but in average can still be described by eq. ^ with 

K 2 cias = o, Kf:: = ln2 . ( l y)2 . > o . (ii) 

7TC44(1 — V) 

It is obvious that such a sample shows 4-th order effective anisotropy with "easy axes" directed along square diagonal 
(45° angle with respect to X axis, as results from condition Kf™ > 0). 

The rotational torque calculated from eqs. JHJ), (J9j) for a typical antifcrromagnet NiO is shown in figs. Q] [2l In 
calculations we use experimental results^ 3 - for a rotational torque around [111] axis taken at room temperature for 
H=4.8 kOe (triangles in fig. [1]). Since the experiment shows strong hysteresis in clock- wise and counter-clock wise 
rotations, the data in fig. [1] are preliminarily averaged over cc-ccw cycle. Theoretical curve (solid line) includes 
some nonzero average torque that may result from inhomogeneity of the sample, and adjusting parameters are x = 
4.35 • 1(T 4 emu/cm 3 , Kf as = 900 erg/cm 3 , Kf™ = 450 erg/cm 3 , a/b = 20. 

Four curves in fig. [^demonstrate possible variation of the rotational torque with the crystal shape. For a large aspect 
ratio (curve (a)) contributions from both 2-nd and 4-th order anisotropy terms are equally important (Jr§ las cx Kf~H), 
the torque curve is composed of sin 29 and sin 46 components. At lower aspect ratio (curves (b), (c)) the 4-th order 
component becomes less pronounced and the amplitude of torque also diminishes. A circular cylinder (a = b) will 
show no shape effect in a single domain state, but a sample with a square cross-section should posses 4-th order 
anisotropy (as seen from curve (d)) regardless of crystalline symmetry. 

Evidently shape-dependence of the magnetic rotational torque for NiO crystal was observed in Ref. [H. The 
authors notice that "a nearly pure sin 46 curve is obtained when the (111) cross section is square". For an arbitrary 
shaped section the experimental curve (fig. [TJ triangles) is satisfactorily fitted with combination of sin 26 and sin 49 
components (solid line). Dashed line in fig. [1] shows theoretical curve with /™^pi anc = 220cos6# erg/cm 3 that could 
be expected in neglecting of shape-induced effects. The presence of sin 49 component in the rotational torque makes 
it possible to exclude the effect of L "freezing" by magnetoelastic strain which may be expected in relatively small 
magnetic field. Additional anisotropy induced by frozen lattice is uniaxial and should be insensible to the variation 
of a crystal shape. 



VI. ANTIFERROMAGNETIC RESONANCE 



The effect of shape-induced anisotropy may also be detected by measuring frequency of the lowest branch of a spin- 
wave spectrum. Antiferromagnetic resonance (AFMR) frequency calculated within a standard Lagrangian technique 
in long-wave approximation is given by the expression 

^afmr = 9\p- (Kg™ + Kg™ cos 26» + 4Kf™ cos 49 - X H 2 cos 2{9 - </>)) , (12) 

where g is gyromagnetic ratio and equilibrium value of 9 may be calculated from Polar diagram of ^afmr("0) 
calculated from (|12p for NiO (g = 2, 5) for a different crystal shape is shown in fig. |U For the aspect ratio close 
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to 1 (dotted line) the shape effect is negligible and magnetoelastic gap in AFMR spectrum is almost isotropic. For 
elongated (dash-dotted and dashed lines) or square-shaped (solid line) samples the AFMR gap should show strong 
2-fold or 4-fold anisotropy. 



VII. SHAPE EFFECT IN ANTIFERROMAGNETIC NANOPILLARS 



Taking into account recent interest to multilayered structures based on antifcrromagnetic materials we consider 
possible shape effects in antiferromagnetic nanopillar that can be a constitutive part of spin-valve structure (see, 
e.gJ^). Typical nanopillar has a form of a very thin (thickness c ~ 3 -r- 10 nm) elliptic cylinder with the pronounced 
in-plane aspect ratio (with a oc 120 nm and b oc 50 nm). 

In this case (c <C a, b) the constants of shape-induced anisotropy may be expressed through the parameter k 2 
1 — b 2 /a 2 that depends on the aspect ratio b/a < 1, namely, 

^das„c [(A') 2 (2-3^) + A.A']J 2 (fc) las ^c 2(A') 2 J 4 (fc) 
2 ~ b' 4c44(l-*/) ' 4an " b ' 3c 44 (l - v) ' 1 j 

where we have introduced the dimcnsionless shape- factors J 2 ,4 as follows 

^ /2 (sin 2 <P + cos2<j)/k 2 )d(j> 



Mk) = f 

Jo 



o v 1 — k sin 



T „. r /2 (1 - 8cos2</>- fc 2 sin 2 ( / 1 + 8cos2(/)/A: 2 )d< ? !) 

■M k ) = / ; = 9 = ■ (I*) 

J a Vl-fc 2 sm <t> 

Dependence of the shape-factors J 2i 4 vs aspect ratio calculated according to eq. (fT¥j) is given in FigJ3] It can be easily 
understood that both shape- induced constants (| 13|) vanish for isotropic sample (b = a). For all the values of aspect 
ratio the 2-nd order term is greater than that of the 4-th order, J 2 > Ji- In the experimentally accessible range of 
values b oc 0.5a J 2 = 0.52 and J4 = 0.3. Thus, characteristic value of the shape-induced anisotropy in the stress-free 
thin film can be of the order magnetoelastic energy Kf^ as oc (A') 2 /C44 multiplied by small factor c/b that varies within 
the range 0.05-^0.3 depending on the film thickness. 

Substantial enhancement of the 2nd order shape-induced anisotropy Kf &s should be expected in the case of mis- 
match between antiferromagnet and substrate lattices. Lattice misfit is a source of rather strong (usually isotropic) 
in-plane stresses = <7™ = cr mf that should be added to the intrinsic stresses J3]). Corresponding (and principal) 
contribution into the effective 2-nd order anisotropy constant takes a form 



Kt-.'f^m. (15 , 

4C44(1 — V 



If u is lattice mismatch and u spon tx A'/c44 is an observable spontaneous strain that occurs in the Neel point, then, 
we can estimate cr mf oc C44W mf (in assumption that the elastic modula of substrate and antiferromagnet are of the 
same order of value), K^ as oc (A') 2 /c 4 4 and hence, 

r ?/ mf 

Kf™ oc Kf} 3 *. 

2 b it s P° n 18 

Substituting typical values of small lattice misfit u mf = 0.005 and large spontaneous striction u s P on = 10 -4 we see 
that even for very thin nanopillars with c/b = 0.05 shape induced anisotropy may be as large as i^2 las = 2.5Kf^ as and 
thus much greater than the "bare" in-plane magnetic anisotropy of antiferromagnet (see table [TJ> . 

Nontrivial relation (fTB")) between the shape of the sample and external stress produced by the substrate may reveal 
itself in switching of the shape-induced direction of easy axis for different substrates. Really, if Kf 1 ^ > 0, then, 
equilibrium orientation of L in monodomain sample is parallel to the ellipse's long axis a (X-direction), as can be seen 
from (J4j) and inset in fig.Q] According to equation (Tl5|) . the sign of i^2 las depends upon the relation between intrinsic 
(A',u spon ) and extrinsic (a ,u ) stresses (or strains). In the case when both substrate and eigen magnetoelastic 
forces of antiferromagnet "work" in the same direction, trying to extend (compress) crystal lattice, product <r A' will 
be positive (in accordance with Le Chatclicr principle) and i^2 las > 0. If we use another substrate which produces 
misfit of opposite sign, i^ 1 ^ < and equilibrium orientation of L will be parallel to the short axis b (^-direction). 

Controlling of spin-orientation by substrate- induced strain was recently observed^ in antiferromagnet CoO. In bulk 
samples CoO is compressed in L direction. When grown on Mg(100) substrate, CoO lattice is expanded in plane and 
experiment shows that Co spin go out of plane. And on contrary, in-plane ordering is observed for Ag(100) substrate, 
which produces slight contraction in the film plane. From our point of view, analogous experiments with nanopillars 
can be very instructive in further study of shape-induced effect in antiferromagnetic crystals. 
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Figure 3: The 2-nd (a) and 4-th order (b) shape-induced anisotropy of thin nanopillar as a function of aspect ratio. Arrow 
indicates the typical value of b/a used in the experiments^. 



VIII. CONCLUSIONS 



In summary, we propose a model that describes an antiferromagnet with the pronounced magnetoelastic coupling. 
The model is based on the assumption that antiferromagnetic ordering is accompanied by appearance of elastic 
dipoles. Due to the long-range nature of elastic forces, the energy of dipole-dipole interaction (destressing energy) in 
the finite-size sample depends on the crystal shape and is proportional to its volume. 

In the "easy-plane" antiferromagnets with degenerated orientation of easy axis the destressing effects may stimulate 
formation of the domain structure and re-distribution of the domains in the presence of external magnetic field. 

The model predicts existence of the shape- induced magnetic anisotropy which corresponds to macroscopic symmetry 
of the sample and can be detected by the magnetic rotational torque and AFMR measurements. 

Crystal shape may be a source of additional magnetic uniaxial anisotropy which produces different effects depending 
of aspect ratio. Below some critical value of a/b shape-induced anisotropy favors formation of the domain structure 
even in the absence of any external field. For large aspect ratio (above the critical value) shape-induced anisotropy 
removes degeneracy of easy-axis in a single domain sample. Energy difference between thus induced easy- and hard- 
directions depends upon a/b. Square cross section sample should acquire 4th order anisotropy (irrespective to the 
crystal symmetry in easy plane) . This opens a possibility to control macroscopic properties of the sample varying its 
shape. 

The shape of antiferromagnetic nanopillars imbedded into structures with lattice mismatch may be a principal 
source of magnetic anisotropy. This fact should be taken into account in engineering of spin-valve devices. 
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Figure 4: Angular dependence of AFMR frequency vs magnetic field orientation for samples with different aspect ratio: a/b = 20 
(dash-dotted line); 3 (dashed line); 1.2 (dotted line). Solid line corresponds to a square-section sample. 
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